Approaching finite-temperature phase diagrams of strongly correlated materials: 

a case study for V2O3 
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Examining phase stabilities and phase eqnilibria in strongly correlated materials asks for a next 
level in the many-body extensions to the local-density approximation (LDA) beyond mainly spec- 
troscopic assessments. Here we put the charge self-consistent LDA+dynamical mean-field theory 
(DMFT) methodology based on projected local orbitals for the LDA-I-DMFT interface and a tai- 
lored pseudopotential framework into action in order to address such thermodynamics of realistic 
strongly correlated systems. Namely a case study for the electronic phase diagram of the well-known 
prototype Mott-phenomena system V2O3 at higher temperatures is presented. We are able to de- 
scribe the first-order metal-to-insulator transitions with negative pressure and temperature from the 
self-consistent computation of the correlated total energy in line with experimental findings. 

PACS numbers: 71.45.Gm, 71.45.Lr, 71.30.+h, 73.20.Mf, 71.15.Mb 



I. INTRODUCTION 



The first-principles computation of phase diagrams at 
finite temperature T for multi-component materials sys- 
tems is a quite formidable challenge. Although there are 
very successful (semi-) empirical methodologies to com- 
pute the thermodynamics of binary (or higher) realistic 
systems, most notably the CALPHAD approach^i the ca- 
pability of predicting phase diagrams by starting from an 
ab-initio quantum-mechanical level has a rather strong 
appeal to many theorists. Dating back to the pioneer- 
ing work in this research area by Hume-Rothery in the 
1930's through empirical rules based on atomic sizes and 
electronegativities;^ the field has reached quite a level 
of sophistication. After the extension of Hume-Rothery's 
original ideas by Miedema and coviforkers^ via addition- 
ally introducing the electronic charge density in the de- 
termination of the formation energy, approaches build on 
density functional theory (DFT) in the Kohn-Sham rep- 
resentation^ eventually have taken over and have been 
dominating the research on atomistic phase-diagram cal- 
culations^"— since the mid 1980s (see e.g. Ref. |l3 for a 
review) . 

However, materials systems with less-screened 
Coulomb interactions among the electrons of the order 
of or larger than the bandwidth W have remained so far 
out of reach. Conventional representations of DFT, via 
e.g. the local-density approximation (LDA)^ are not 
capable of accounting for the effects of strong electronic 
correlations. Phase transformations at finite T either of 
pure electronic kind or driven by electronic correlations 
are usually not dcscribablc solely within Bloch band 
theory. On the other hand, many novel materials 
which are technologically promising, because of e.g. 
enhanced response behavior, display signatures of strong 
correlations. Furthermore, even well-known allotropes 
of transition metals or prominent transition-metal 
alloys with or close to magnetic order (like e.g. the 
iron-aluminum systemii) are rather difficult to model 
within standard LDA(-likc) approaches due to the lack 



of explicit many-body correlation effects. 

The combination of LDA with the dynamical mean- 
field theory (DMFT), the so-called LDA-^DMFT ap- 
proach, nowadays prosperously allows to include the ef- 
fects of strong Coulomb interactions in realistic solids 
(see e.g. Ref. [ij for a review). Note that one may easily 
use the generalized-gradient approximation (GGA) for 
the DFT part, but as far as it concerns explicit strong 
correlation effects the difference to the LDA approach 
are usually negligible. Yet there arc only few imple- 
mentations that handle the LDA-I-DMFT formalism in 
a charge self-consistent framework, i.e. accounting for 
the feedback of the local electronic self-energy onto the 
charge density that determines the Kohn-Sham effective 
potential, until self-consistency of the complete interact- 
ing charge density is achieved J^"— Additionally, in order 
to account for competing strongly correlated phases at 
elevated temperatures there are high demands on the ac- 
curacy and generality of the underlying band-structure 
methodology as well as the utilized DMFT impurity 
solver. 

As to explicit realistic phase-competition studies, there 
is prominent work within that scope mainly in the area 
of /-electron compoundsJ^i^ In a recent LDA-I-DMFT 
study, Leonov et al. provided a quite successful modeling 
of the high-temperature bcc-to-fcc transition in iror^ us- 
ing the powerful quantum Monte Carlo method for the 
impurity solution. However, there the charge selfconsis- 
tency was neglected and total energies have been calcu- 
lated in a post-processing scheme. 

In the present work we want to review the current 
state of the art formalism for handling the charge self- 
consistent LDA-I-DMFT method with the direct calcula- 
tion of the correlated total energy. The full approach 
is applied to the key features of the electronic phase 
diagram of the famous sesquioxide V2O3 (see Fig. [T]) 
above room temperature. When revealing the phase 
boundaries between metallic and insulating phases with 
T and pressure p, we neglect the effect of chemical dis- 
order as well as explicit electronic entropy contributions. 
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FIG. 1. Experimental phase diagram of V2O3 from Refs. Il£ 
and [23. 



The former is not expected to play a vital role in the 
present study since we are interested in the stoichiomet- 
ric system well below possible disordering/melting tem- 
peratures. Extending the strongly correlated formalism 
to composition-dependent phase diagrams is in principle 
possible, e.g. via the so-called stat-DMFT methodi^ As 
it will be shown, the present results in the more restricted 
scheme arc already encouraging and it should be mainly 
a matter of time when such new and further elaborate 
techniques will become a generic tool in the context of 
first-principles phase-diagram evaluations. 



II. THEORETICAL FRAMEWORK 

The charge-self-consistent scheme of the LDA-f DMFT 
framework is technically and computationally rather de- 
manding and only few implementations thereof exist up 
to nowfii^i^ For the DFT(LDA) part of the calculations, 
we here employ a mixed-basis pseudopotential^i (MBPP) 
as well as a projector-augmented wave^ (PAW) imple- 
mentation. 

The explicit many-body effects are treated within 
a multi-orbital Hubbard model including also all non 
density-density local interaction terms, i.e. we use the 
complete Coulomb matrix according to full rotational 
symmetry. This amounts to a parametrization with two 
interaction integrals, namely the Hubbard U and the 
Hund's exchange J. The interacting problem is solved 
within DMFT using a continuous-time quantum Monte- 
Carlo (CT-QMC) impurity solver in the hybridization- 



expansion formulation^^ as implemented in the TRIQS 
package 1^ 

In short, the chargc-sclf-consistcncy condition can on 
equal footing be understood as a way of improving the 
density functional of DFT to include more explicit cor- 
relation effects as well as finding a realistic and consis- 
tent effective single-particle part of DMFT via combining 
both formalisms to a new, complete cycle, which is sum- 
marized in the following sections. 

Before that, however, it is vital to note that on gen- 
eral grounds this LDA-f-DMFT formalism is manifestly 
temperature dependent and therefore in principle ideally 
suited for thermodynamic problems. This is in stark con- 
trast to the standard extension of the Kohn-Sham formal- 
ism towards finite T via the Mermin theory of including 
the proper Fermi function for the electronic statesi^ Here 
the full impact of temperature on the many-body level, 
including e.g. the effective disappearance of Bloch-like 
quasiparticle states and thus the localization due to the 
loss of coherency at large T, is properly taken care of. 



A. Projected local orbitals 

The DFT(LDA) method utilizes an orbital- 
independent representation of the effective-single 
particle Hamiltonian for the electronic structure result- 
ing in Bloch-Kohn-Sham (KS) wave functions for the 
solid-state electronic structure. On the other hand the 
DMFT equations make use of a local correlated subspace 
in order to include the effects of strong Coulomb cor- 
relations in condensed matter. Thus the first step that 
has to be peformed when interfacing the DFT(LDA) 
and the DMFT technique is the extraction of a suitable 
correlated subspace C starting from the complete Hilbert 
space of Bloch KS band states. This is done in the 
projected-local-orbitals (PLO) schemei^"— For com- 
pleteness we here briefly review the methodology in the 
context of charge-self-consistency and the total-energy 
calculation. More details may be found in Ref. [Sll 

Normalized orthogonal projections onto chosen local 
orbitals with character m and centered at site R may be 
defined via 
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where the \^kv) for wave vector k and band v are chosen 
to be a subset W of the Bloch states of the original LDA 
treatment and O describes the overlap matrix written as 






l*fcJ(*feJXfl'. 



(2) 



The set of states {|x_Rm)} together with the energy win- 
dow W define the correlated subspace C , chosen such that 
the problem of strong local Coulomb interactions is ad- 
equately represented. Conveniently, C is adapted to an 



available localized basis used in a given band-structure 
code, i. e. linear combinations of the mixed basis within 
the MBPP framework or the partial waves of PAW. 

Using the projections ([1]), one can construct the one- 
particle Green's function within the truncated Bloch 
space W via the double-counting corrected local DMFT 
self-energy (written in Matsubara frequencies uj„={2n + 
1)77/03^), which is assumed block diagonal in the corre- 
lated sites, written as 



A^^'(jc.„) = (i™P'^(*c.„) -^'^^^^)<5w 



(3) 



This Bloch Green's function is thereby connected to the 
correlated subspace through an upfolding procedure, i.e.. 
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In this equation, cr denotes the diagonal matrix of 

Kohn-Sham eigenvalues for the Bloch states and fi marks 
the chemical potential (see section Hi D I for details). This 
Green's function can then be downfolded to the corre- 
lated subspace, enforcing the DMFT self-consistency con- 
dition (proper normalization of the fc-sum here and in the 
following is understood) 

g-P'^(zc.„) = 5]£-(fc) • g'\lUJn,k) ■ ^\k) . (5) 



This impurity Green's function can then be used to sup- 
ply a new DMFT bath Green's function 

(£^^)~'(*w„) = (£™P'^)-^(zw„) + g"^P^^(*w„) (6) 

that enters the impurity solver yielding eventually an up- 
dated impurity self-energy until convergence is achieved. 
The outlined iterative scheme marks the usual DMFT cy- 
cle without charge-self-consistency, since it then works as 
post-processing scheme to a once computed set of Kohn- 
Sham objects {e^s, ^^^ J. 



B. Expressing charge densities 

The fundamental step in the self- consistent combina- 
tion of DFT(LDA) and DMFT is provided by the ex- 
pression of the basic quantities of each method in terms 
of the basic quantities of the other method. Namely, 
charge density for DFT and one-particle Green's func- 
tion for DMFT. For this purpose, we define a Kohn-Sham 
Green's function through 



G (ia;„,fc) = (iw„ + ^Ks)l 



,KS 



(7) 



Note that this function is in general different from the 
DMFT "Weiss-Field" ^o- The choice of /iks is described 



in detail in section Hi Dl In the following the band indices 
vv' live in the truncated Bloch Hilbert space W and we 
drop for convenience the site index R. Generalization 
of the formulae including the latter is straightforward 
and since the charge density is additive, contributions 
from (supposingly weakly correlated) bands outside W 
are most easily taken care of. The trace of G (iw„,fc) 
expressed in the Bloch basis is nothing else than the 
charge density of a standalone KS-LDA calculation which 
reads 

P^'d) = ^ E (^l*^'') G^i^^n.k) {^,.,\r) , (8) 

knuu' 

with p^l/k^T as the inverse temperature. A very simi- 
lar form can be found for the charge density from a post- 
processing DMFT calculation, i.e. 



knuf' 

Thus the difference p'{r) = piz) — P^^{l) is given by 
1 
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As described in Ref. [SSJ, the object AN (k) can be rewrit- 
ten as 

Aiv(fc) = ^ E{^'^^(*^"'^) ■ (£^® • ^i(*^») ■ £® 

n 

-(/^-MKs)i)-a"(*^n':^)}- (11) 

Therewith a simple representation of the total charge 
density including self-energy effects beyond LDA is pro- 
vided, reading 

k^vv' 

(12) 
and to be used and manipulated in a given DFT-based 
band-structure code. Here /(e) denotes the Fermi- 
distribution function and lkv=^kv—PKS- Hence the inclu- 
sion of the DMFT self-energy renders it necessary to not 
only incorporate modifying terms diagonal in the Bloch 
states, but also off-diagonal contributions. The problem 
of truncating the whole Bloch space to a subspace W 
therefore reduces to taking into account the correct set 
of bands in each summand. Details on the implemen- 
tation thereof in the different KS basis sets are given in 
appendix [Xj 



C. Self-consistency condition 



E. Total energy 



The aim of charge-self-consistency is to include DMFT 
self-energy effects in the charge density, so that p^^(r) 
and p{r) can obviously not be the same quantity. Instead, 
it is most instructive to use the (spectral density-) func- 
tional approach by Savrasov and Kotliar,— incorporating 
both one-particle Greens 's function and charge density. 
Extremization thereof with respect to the charge density 
p{r) basically yields the Kohn-Sham equations in which 
the correlated charge density p(r) is used as an input for 
the effective potential V^ff [p(zi)]i i-e. 



f + Vos[pir)]-ek, \-9k.) =0 



(13) 



Similarly, extremization with respect to the Green's func- 
tion yields the usual expression for the DMFT self-energy. 
Thus the complete cycle can be constructed as follows: 
From an initial DFT(LDA) calculation, perform conven- 
tional DMFT steps and compute a correlated charge den- 
sity p{r) as given by eq. (fT2|) . That charge density is 
reinserted into the band-structure code (for this step, 
knowing the elements of AN (k) is sufScient) and new 
effective single-particle wave functions |^fci/) are com- 
puted using eq. (fT3|) . These finally enter eq. ([T]) to build 
a new correlated subspace for DMFT (keeping the set 
(IXijm)} unaltered). This enlarged cycle is iterated until 
full charge-self-consistency is reached, i. e. charge den- 
sity and self-energies (and thus the matrix AiV'-' ) remain 
constant with iterations. Our experience with the cycle 
shows that convergence is robust, a linear mixing scheme 
is often sufficient. 



D. Chemical potential 

As usual, the chemical potential p is adjusted such 
that the resulting total charge density ptot(z:) holds the 
correct total number of electrons TVtot and is enforced via 



dr ptotir) = iV^^w + ^ X] "^^ (^'^n, fc) 



m. 



(14) 



In Ref. |33| it is argued that this charge neutrality condi- 
tion is imposed on ptot(i:) only (not on p^^{r)). Thus the 
parameter /xkSj which is deffiied to be the energy up to 
which the Kohn-Sham states of the DFT(LDA) part are 
filled (with the proper Fermi-distribution function), can 
be chosen to be equal to the chemical potential p, which 
would mean that the integrated charge density changes 
due to correlations in eq. (|12p vanishes. 

However, in order to clarify relations between some 
of the quantities that occur in the formalism, it can be 
useful to choose ^ks such that the DFT(LDA) part of 
the calculation is already charge neutral. It can easily 
be shown that with the correction term of eq. (|lip , this 
choice does not affect p{r). Anyway, it is important to 
note that /iks itself has no physical interpretation in the 
enlarged LDA+DMFT framework. 



In order to obtain total energies within the charge-self- 
consistent formalism, the spectral density functional ap- 
proach is applied again. Based thereon, the total energy 
may be computed froroi^ 



E, 



LDA+DMFT 



= E 



LDA 



E'l'^N^^ 



(Hn) - E, 



dc 



(15) 
Note that this expression relies on the fact that the 
adapted diagonal basis for the Kohn-Sham single-particle 
Hamiltonian is utilized, which yields eigenvalues e^^j. 
Several approaches are possible to obtain the expecta- 
tion value of the two-particle Hamiltonian {Hu). Here 
we choose to apply the Galitskii-Migdal formula;^ As 
shown by Bochnke et al^, the quality of the numerical 
data can be improved by choosing a suitable basis set, 
i. e. Legendre polynomials, for the represenation of the 
one-particle Green's function. Note that although this 
LDA-I-DMFT total energy is temperature-dependent, we 
here keep the 'energy' notion, since for a well-defined 'free 
energy' a clear definition of an entropic part would be in 
order. 



III. THE V2O3 SYSTEM 

The vanadium sesquioxide V2O3 belongs to the most 
prominent strongly correlated compounds and has al- 
ready been subject to many theoretical effortsi^"— At el- 
evated temperatures it orders in the corundum structure 
in which there are V-V pairs along the crystallographic 
c-axis and a honeycomb lattice appears in the a;y-plane 
(see Fig. [2]) . The V ions reside within an octahedron of 
oxygen ions, respectively, building up a trigonal crystal 
field for the transition-metal ion. Thus the low-energy 
t2g orbitals of the V(3d) shell are split into an oig and 
two degenerate e' orbitals. Formally the vanadium ion 
has the 3(P valence configuration, i.e. is in the V'^^ oxi- 
dation state. The t2g orbital degrees of freedom appear 
to play a central role for the intriguing physics of this 
transition-metal oxide and Castellani et al^ were the 
first to provide a detailed account of the complex corre- 
lated electronic structure in these local terms. 

The finite-temperature phase diagram, taken from the 
original work of McWhan and coworkeraiSi^^ is shown 
in Fig. [1] In this work, V substitution by Ti (Cr) im- 
plies positive (negative) pressure. It displays three major 
phases, namely paramagnetic metallic (PMM), paramag- 
netic insulating (PMI) as well as an antiferromagnetic 
insulating (AFI) regime. The transition to the latter 
AFI phase at lower temperatures is also associated with a 
structural transition to a monoclinic low-symmetry struc- 
ture.— Furthermore, additional phase-diagram studies in 
the vanadium-deficient V2-j,03 regime revealed the ex- 
istence of a metallic spin-density-wave (SDW) phase4^ 
Interestingly, that phase appears not to be a form of 



Mi. 







FIG. 2. (Color online) Corundum structure of V2O3 viewed 
along the z-axis (left) and along the j/-axis (right), with de- 
picted V ions (large grey) and O ions (small red/dark). 



precursor to the much more extended AFM ordering of 
the insulator. On the contrary it seems that the mag- 
netic short-range order within PMM and PMI is closer 
to the SDW ordering. Therefore the magnetic ordering 
in the AFI phase may be closely related to the structural 
change, even involving additional orbital ordering.— Due 
to this additional complexity in connection with the AFI 
phase we concentrate in this work only on the phase equi- 
librium between PMM and PMI. A complete description 
of the V2O3 phase diagram including the magnetically 
ordered phases wiU be postponed to future studies. 



A. LDA characterization and local projections 

A thorough first-principles DFT(LDA) description of 
metallic V2O3 at normal pressure and without doping has 
initially been given by Mattheissi^ Here we only summa- 
rize the most relevant features as they evolve from our 
MBPP investigation. Fig. [3] shows the LDA band struc- 
ture along high-symmetry lines within the first Brillouin 
zone. It is evident that the bands group in a way canon- 
ical for many transition-metal oxides. The larger block 
below the Fermi level ep in the range [—8, —4] eV is dom- 
inated by oxygen 2p orbital weight, while the unoccupied 
block within [2,4] eV stems major ly from vanadium eg 
orbitals. This encoding is visualized in the dcnsity-of- 
states (DOS) plot of Fig.|3]from local projections onto the 
symmetry-adapted cubic-harmonic angular-momentum 
channel of the V{3d) basis. The band manifold of width 
W^2.6 eV around ep is mostly composed of aig and e'g 
orbitals with only minor inter-mixing of V(eg) and 0(2p). 
Note that especially the aig character shows a promi- 
nent bonding-antibonding signature in the DOS of this 
low-energy region. 

From the LDA result the set of local orbitals (IXi^m)} 
to be utilized in the local projections defined in eq. (P) 
are here chosen to be given by the linear combinations of 
pseudized atomic V(3d) functions that diagonalize the 



orbital density matrix on each of the four symmetry- 
equivalent vanadium ions within the unit cell. This is 
often referred to as the crystal-field basis, whereby here 
the aig level is higher in energy than the e' one. Figure^] 
exhibits the local DOS obtained from the projections us- 
ing for W the low-energy t2g manifold. In the following 
we will only concentrate on these minimal projected lo- 
cal orbitals and will not elaborate on the possible cases 
of larger energy windows W, i.e. such ones that also 
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FIG. 3. LDA band structure of V2O3. 
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FIG. 4. (Color online) Total LDA density of states of V2O3 
and local DOS within the symmetry-adapted V(3d) basis with 
range rc=2.0 a.u. 




FIG. 5. (Color online) aig and e'g LDA-DOS on the basis 
of the orthonormalized projected local orbitals, whereby the 
range of the {|Xflm)} was also rc=2.0 a.u. 





1 



9 



-n.i 


W 


^ 


^ 


1 


-i).;5 


% 









-{).-! 


h 


A 


S 


1 



(Color online) LDA bond charge densities 



^atomic ^^^ for V2O3 within the a;z-plane (left) and 



FIG. 6. 
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the xy-plane (right). Note that in the latter case the two cen- 
tral V ions are not at the same height and therefore appear 
different. 



cover the high-energy occupied/unoccupied band mani- 
folds. Note that the local orbital DOS in Fig.Hand Fig.E] 
differ on principle grounds (see also appendix lA3[) . In the 
latter case it is computed in the basis of orthonormal- 
ized orbitals, whereas in the former case it is calculated 
from projections onto angular-momentum channels with- 
out proper final normalization (i.e. no radial orbital func- 
tion involved). From the orthonormalized projected lo- 
cal orbitals one retrieves the occupations nLDA(aig)=0.57 
and nLDA(eg)=1.43 (summed over both e'g orbitals), re- 
spectively. These values differ by about 0.05 electrons 
towards stronger orbital polarization compared to the 
numbers presented in Ref. |4J from a Wannier construc- 
tion within the A^th-order mufhn-tin-orbital method4i 

The stronger ionic character of the transition-metal 
oxide compared to ordinary metals or intermetallic com- 
pounds becomes clear from the plot of the bonding charge 
densities in Fig. [B) The latter function is defined as 
the difference between the crystal valence charge density 
and the superposed atomic valence charge densities. The 
charge transfer from vanadium to oxygen is obvious, but 
also the expected charge accumulation in the interstitial 
region is visible. 



B. LDA+DMFT: finite-temperature phase 
equilibria 



The LDA-only description does not account for a 
metal-insulator transition (MIT) in V2O3. We model 
the interacting problem on the many-body level within 
a three-orbital (aig, 2e' ) multi-site (four V ions in the 
corundum unit cell) generalized Hubbard model employ- 
ing the complete rotational invariant Coulomb interac- 
tions on the local level. For the parametrization of 
the Coulomb integrals we choose the values C/=5 eV 
and J=0.93 eV, as already utilized in earlier simpli- 
fied LDA-hDMFT studies for VaOs.i^ The following re- 
sults are obtained from our MBPF-code interface of 
LDA+DMFT. A second implementation within the PAW 
approach is also briefly discussed in the appendix. 



1. Lattice expansion and temperature variation for fixed 
c/a ratio 

The MIT between the PMM and the PMI phase with 
negative pressure is depicted in Fig. [7] for selected T. 
While the p<0 scenario is realized experimentally via Cr 
doping, it is here provided in a simple way by increasing 
the lattice constant a starting from its experimental^ 
equilibrium value ao=4.95A. Notably we also first keep 
the c/a ratio fixed to its value c/a=2.83 at ambient pres- 
sure and temperature. The effect of relaxing that ratio 
will be discussed in section fill B 31 For the moment, this 
approximate theoretical approach proves to be sufficient 
to describe the key features of the V2O3 phase diagram 
above room temperature. 

One could connect this approach to the physical pres- 
sure p in a simple way by defining p=—dE/dV. It is 
seen that the theoretical formalism reveals the pressure- 
induced first-order MIT with the correct positive sign of 
the slope OTmit / dp from experiment (compare Fig. [Ij. 
However the changes that occur in the lattice constant 
with T along the phase boundary are non-surprinsingly 
larger (roughly by a factor 5-10) than in experiments^ 
The neglect of electronic and phononic entropy contribu- 
tions (and presumably also non-local correlations) may 
be blamed. Nevertheless the change of curvature of the 
respective total energies elucidates the expected soften- 
ing of the lattice with increasing temperature from the 
decrease of the bulk modulus B^^y^. Fig. [8] displays the 
tie-line construction for the first-order transition between 
the metallic and insulating phase. The volume jump at 
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FIG. 7. (Color online) MIT with negative pressure, i.e. in- 
crease of the lattice constant for various temperatures within 
LDA+DMFT. Left: total energy -Elda+dmft normalized to 
the value at the equilibrium lattice constant, right: energies 
with respect to the LDA energy -Elda for each given lattice 
constant. 
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FIG. 8. (Color online) Tie-line construction for the first-order 
MIT with negative p for r=387 K, 290 K, 232 K (from top 
to bottom). 



the transition becomes obvious, albeit again the retrieved 
pressure from the slope of the common tangent is about 
an order of magnitude too high compared to experiment. 
But note that the increased stability range of the PMI 
phase with volume at larger T is a direct result of the 
calculations. For rcrit^400 K the line of first-order tran- 
sitions at negative pressure exibits a solid-solid critical 
end point and a continuous path from the metallic to the 
insulating regime opens. One may already recognize that 
the shifted quasiparticle peak at the lower gap edge in the 
spectral function of the insulator has vanished for T=387 
K, which may signal the immediate strongly incoherent 
regime close to the critical end point. 

The pressure-dependent investigation describes more- 
over the increase of correlation strength when approach- 
ing the critical p deep from the metallic regime. This 
is documented in Fig. [5] where we plot the local spec- 
tral function with increasing lattice constant. Strong 
transfer of spectral weight from the low-energy region to 
the high-energy Hubbard bands is observable and marks 
the evolution towards the MIT with negative pressure at 
constant interaction strength. Furthermore a shifting of 
the dominant quasiparticle peak towards the Fermi level 
with increasing the lattice constant may be recognized. 
Of course, the growing lattice distances also weakens the 
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FIG. 9. (Color online) Evolution of the metallic spectral func- 
tion with negative pressure for r=387 K. The inset shows a 
blow-up of the features around the Fermi level. 



metallic screening and therefore should lead on simple 
grounds to an effective increase of the mutual Coulomb 
interaction between the electrons. Such an effect is here 
at least in simplest terms describable by the charge-self- 
consistent reaction to the applied lattice expansion. Note 
however in this context that recent photoemission studies 
point to a rather constant U value on the different sides 
of the phasc-boundaries;^ 

In order to trace the MIT with temperature in 
some detail, we plot in Fig. [10] the intrinsically T- 
dependent LDA+DMFT energy now at fixed elongated 
lattice constant. Starting from the low-temperature 
metallic regime, £'LDA+DMFT(r) deviates from a sim- 
ple functional behavior at T~^270 K, displaying an 
overall double-parabolic structure. The MIT takes 
place around Tmit^SIO K in good accordance with 
the expected experimental region of the temperature- 
induced PMM— PMI transition. Surprisingly, the energy 
^lda-i-dmft(2^) of the PMI phase appears rather flat in 
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FIG. 10. (Color online) MIT with temperature for the fixed 
lattice constant a=1.08ao, i.e. at finite negative pressure. 



the temperature regime [350, 450] K, which is just in the 
neigborhood of the experimental critical end point. 



2. Charge densities and orbital resolution for fixed c/a ratio 



f> ^f *> i^ 




So far we concentrated on the integral impact of the 
electronic correlation on the finite-temperature prop- 
erties of the V2O3 system. However, for a deeper 
understanding of the underlying driving forces it is 
also important to shed light on the possibly distinct 
behavior of individual microscopic degrees of freedom 
and most notably on those of orbital kind. Concern- 
ing the distinct orbital occupations with temperature 
and negative pressure no dramatic effects occur in the 
correlated electronic structure. In line with previous 
postprocessing studies^^i^ within charge-self-consistent 
LDA-I-DMFT the aig orbital filling is generally reduced 
compared to the LDA value (and correspondingly the 
e' filling is increased). For the equilibrium volume and 
r=232 K the numbers write as nDMFT(Q'ig)=0.48 and 
?^DMFT(e'„)=1.52. A real-space discrimination of these 
orbital fiUing differences between LDA+DMFT and LDA 
on the basis of the respective charge densities is displayed 
in Fig.[TT]for T close to the critical end point. In general, 
a localization effect takes place, whereby charge from the 
interstitial region is transfered closer to the atomic sites. 
Thereby the aig orbital (pointing roughly along the z- 
axis) looses charge in the correlated electronic structure, 
whereas the e' orbitals gain4^i^ In the insulating regime 
less charge is transfered from the interstitial part, but 
note that here the lattice constant is also larger. 

Within the correlated scheme, increasing the negative 
p, i.e. enhancing the lattice constant, yields a slight filling 
increase of the aig orbital. When raising T for a=1.08ao, 
so that the system shows a temperature-induced MIT, 
the same trend occurs. The same effect with tempera- 
ture was already theoretically observed in Ref. |4J. Thus 
the calculation reveals an increase (decrease) in the oc- 
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FIG. 12. (Color online) Difference pSSFT(r)-PDMPT(r) 
between the LDA-I-DMFT charge densities in the metallic 
(r=232 K) and in the insulating (T=387 K) regime for 
a=1.08ao within the sz-plane (left) and the xy-p\a,ne (right). 
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FIG. 13. (Color online) Orbital-resolved local spin suscepti- 
bilities with temperature (left) and negative pressure (right). 



cupation (of the order of a few percent) for the aig (e' ) 
orbital in the insulator compared to the metal. It is 
again instructive to visualize directly the changes in the 
self-consistent correlated charge density. Figure [12] de- 
picts the differences in the LDA-f-DMFT charge density 
Pdmft(zi) at the different temperatures associated with 
the metallic and the insulating phase. It is seen that 
here the charge transfers are marginal, mainly showing 
the aig orbitals gaining some charge against the e' with 
T. As an effect of the temperature raise, the interstitial 
inbetween the V ions appears to loose some small weight, 
in the spirit of effective localization at high T. 

Finally, we turn to a brief look on the magnetic re- 
sponse. The orbital-resolved local spin susceptibility Xj 
plotted in Fig. [T31 shows for both orbital contributions 
the expected Curie-Weiss-like tail at higher tempera- 
tures, but for fixed lattice constant a non-monotonic be- 
havior below r~270 K emerges, i.e. in the same range 
where the non-trivial characteristic in £'LDA+DMFT(r) 
was observed. The quenching of Xai in that regime needs 
further study and might be interesting in the context of 
orbital-ordering in the low-temperature AFI phase. As 
expected, for constant T the increase of the lattice con- 
stant leads also to nearly constant x in the PMI phase. 



3. Effect of relaxing the c/a ratio 



FIG. 11. (Color online) Differences of the electronic charge 
densities Pdmft(i:)— Plda(ii) at T=387 K for the metal with 
a=ao (left) and the insulator with a=l.lao (right) within the 
a;2-plane (top) and the xy-plane (bottom). 



From experiment it is known that the c/a ratio is low- 
ered when passing from the metal to the insulator with 
negative pressure i^i^i^ In order to account for that ef- 
fect we relaxed c/a for each volume V within charge-self- 
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FIG. 14. (Color online) Left panel: Total enery vs. volume 
for fixed c/a (solid) and relaxed c/a (dashed). Right panel: 
Relaxed c/a values vs. volume (dotted). Black lines mark the 
metallic solution, red (grey) lines the insulating one. All data 
for T=387 K. 



consistent LDA+DMFT at T=:387 K by computing the 
total energy for selected c/a values and finding the min- 
imum E{V, c/a) via a polynomial fit to the data points. 
Note that c/a also varies substantially with tempera- 
ture,— however as a proof of principles we here only 
followed its evolution with expanding unit-cell volume. 
The results are shown in Fig. [M] No dramatic effect 
results in the global energetics, however the respective 
energy gain, especially in the PMI phase, is clearly vis- 
ible. While in experiment the c/a value varies within 
the interval [2.78,2.88],-^ the given range [2.70,2.80] is 
somewhat larger from the calculations but still within 
the right ballpark. In accordance with the experimen- 
tal data the PMM phase has a larger c/a than the PMI 
phase. Interestingly, the minimum ratio is just reached 
in the transition region with respect to the volume, i.e. 
c/a{V) becomes most soft close to the negative-pressure 
driven MIT. Of course, at the latter first-order transition 
a jump in c/a will take place in line with the volume 
jump corresponding to the tie-line construction. 

A relaxation of c/a also affects the trigonal crystal- 
field splitting At between aig and e'g orbitals. Because 
of the lowering of c/a with negative p and hence a weak- 
ening of the distortion of the VOg octahedra a reduc- 
tion of At is expected. Indeed the calculations reveal a 
shrinking of Af from the PMM to the PMI phase. Note 
that thereby At is determined from the sole Hamiltonian 
part of the interacting problem stemming from the local- 
orbital projections acting on the LDA-t-DMFT converged 
Bloch states. A second contribution to the resulting ef- 
fective crystal-field splitting Aoff is given by the real part 
of the self-energy difference between aig and e'g at zero 
frequency, i.e. Acff=At-h3fiEaig(0)-3?Ee^(0) .ii However 
the latter difference differs only little between PMM and 
PMI and thus no dramatic changes in the occupations 
occur. 



Contrary to that, former post-processing LDA+DMFT 
studies with an interaction-driven picturing of the differ- 
ence between PMM and PMI phase resulted in rather 
strong orbital-polarized solutions for the insulating case. 
That went along with an enhanced value for Aotr due 
to large self-energy effects. Note that in the charge-self- 
consistcnt LDA-I-DMFT framework a clear-cut separa- 
tion into one-body and many-body contributions to Acg 
is not that simple anymore. The reason is that during the 
self-consistency cyle the hoppings also change because 
of the self-energy effects, contrary to post-processing 
LDA-f-DMFT. In our work we did not change the in- 
teraction parameters within the different phases, in line 
with recent experimental worki^ We however tried to 
change our global ([/, J) parameters in a certain range 
to look for the possibility of orbital-polarized solutions. 
Because of the complexity of the problem depending on 
the interplay of hopping and many-body effects within 
an evoluting crystal-structure evolution temperature and 
pressure, we may not exclude a certain setup that allows 
for orbitally polarized solutions. Yet within our studies 
we did not find clear evidence for this behavior. 



IV. DISCUSSION AND SUMMARY 

In this work the advancement of the LDA+DMFT 
methodology, namely the implementation of a complete 
charge-self-consistent scheme with total-energy calcula- 
tion, build on a pscudopotential band-structure code 
using plane waves and localized functions, was docu- 
mented and the principle formalism of interfacing LDA 
and DMFT utilizing projected local orbitals was re- 
viewed. In the calculations we observed that the charge- 
self-consistent framework, at least when using a minimal 
energy window for the projected local orbitals, leads for 
fixed values of the Hubbard U to somewhat smaller elec- 
tronic correlations than the elder post-processing scheme. 
Also orbital polarizations are generally slightly weaker in 
the new complete methodology. This outcome might be 
not that surprising since the now possible reaction of the 
charge density to the self-energy effects may lead to ad- 
ditional screening effects. Further more detailed investi- 
gations of the charge-self-consistent technique are needed 
in this respect and shall be hereby stimulated. 

As a proof of principles, this approach renders it pos- 
sible to describe the first-order character of the MIT in 
the challenging V2O3 system induced by negative pres- 
sure and temperature in accordance with the experimen- 
tal phase diagram. The methodology is in the position 
to describe the PMM/PMI phase boundary in a quali- 
tative correct manner. Yet the quantitative agreement 
concerning structural data and pressure in the transition 
region is still not perfect. The neglect of electronic and 
vibrational entropy terms may be a probable cause, also 
since the absolute value of the critical negative pressure 
|p|<lGPa is rather small compared to other pressure- 
driven MITs^ 
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At first sight, a qualitative difFcrcncc concerns the re- 
spective orbital {aig, e' ) fillings in the PMM and PMI 
phase. Whereas polarized x-ray absorption measure- 
ments together with multiplet calculations point to an 
increased orbital polarization towards less filled aig in 
the insulating regime,— our calculations result in a ten- 
dency towards slight orbital balancing in the PMI phase. 
The named stronger orbital polarization in the insulating 
regime is indeed verified in post-processing LDA-I-DMFT 
for fixed lattice structure, larger U and constant tM^ 
It results there from the increased effective crystal-field 
splitting due to the strong electronic correlations. How- 
ever it is important to note, that in the experimental 
part of the work by Park et al— the negative pressure 
regime was realized by non-isovalent Cr substitution for 
V. Recent work by Rodolakis ct aij^— compared that 
Cr doping-driven (and thus implictly negative-pressure 
driven) MIT with a true pressure-driven MIT. The lat- 
ter scenario was realized by increasing pressure on insu- 
lating (Vo.92Cro. 28)203. It was observed that with true 
applied pressure the orbital occupations hardly vary at 
the MIT, in good accordance with our results. Hence 
rather strong orbital polarization in the PMI phase ap- 
pears to be bound to the doping-driven realization of the 
insulating phase. Further theoretical studies that explic- 
itly treat the chemical doping within an e.g. supercell 
approach, should clarify this issue. The former post- 
processing LDA-I-DMFT treatments followed the route 
of mainly interaction-driven MIT, while from recent hard 
x-ray photocmission spectroscopy it was concluded that 
the Hubbard U docs not change through the MITi^ 

Of course when it comes to the question of orbital po- 
larization, the competition between crystal-field effects 
and Coulomb correlations depends on the choice of the in- 
teraction parameters (C/, J)^ and on the energy window 
used for the projection onto the correlated subspace. Es- 
pecially concerning the latter, e.g. inclusion of the 0{2p)- 
dominated bands (and/or extending the many-body part 
to a five-orbital sector) will change the notion of the 3d 
orbitals and their overall occupations significantly. More 
detailed studies along those lines, also by utilizing ab- 
initio computed Coulomb integrals, are surely necessary. 

But despite this need the current work renders it clear 
that there is advancement in the LDA-I-DMFT frame- 
work that sharpens the tool for strongly correlated mate- 
rials investigations of finite-temperature phase competi- 
tions on an equal footing with analyses of the involvment 
of local degrees of freedom. 
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Appendix A: Representation of ciiarge densities in 
the Kohn-Sham basis sets 



In this appendix we summarize how to compute the 
given charge density from the matrix AN}-' (as in eq. 
P^ ) in two possible Kohn-Sham basis sets. 



1. Mixed-Basis Pseudopotential (MBPP) 

The mixed-basis pseudopotential (MBPP) ap- 
proach^liS uses normconserving pseuodpotentials^ and 
a combined basis of plane waves and modified localized 
atomic functions (f>~i„^{r) for the representation of the 
crystal wave functions, written as 



*fe^(i:) 



1 



.J2^p^ik+G)r 



Palm't'almV-) 



aim 



(Al) 
where ^g and P^im are the respective expansion coeffi- 
cients for atom a in the unit cell and angular-momentum 
numbers Ivi. The correlated charge density therefore 
consists of three parts p{r)=p'^^\r)+p^'^\r)+p'^'^\'r), cor- 
responding to a plane-wave term, a mixed term and a 
localized-function term. With the abbrevation 



N, 



{k) 



fih.)S..' + ANl^l 



(A2) 



they can be written as follows 



.i{Q.'-G}r 



(A3) 



GG' 



P^'\r) 



p(3)(r) 



J{k+G)r 



^ ^ fci/i/' G 

aim 

kyiy' aim 

a'V m' 



(A4) 



(A5) 



The first term can be evaluated directly by Fourier trans- 
formation of both wave functions to real space. The sec- 
ond and third term, which arc zero in the interstitial 
region (where the localized function have decayed), are 
calculated in a straightforward way in an atom-centered 
basis. 



2. Projector-Augniented Wave (PAW) 

The implementation within the PAW formalism is in 
line with Ref. [13. As shown in Ref. [5^, charge densities 
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FIG. 15. (Color online) The equation of state of V2O3 cal- 
culated with the conventional LDA+U scheme as well as the 
charge-self-consistent scheme using an Hartree-Fock solver for 
the DMFT impurity problem. 



here break down also into three parts, namely a plane- 
wave part p{r) and a one-centre term from partial waves 
Pr{l) and from pseudo- partial waves p}j(r) per atom at 
site R (omitting core densities, which arc not affected by 
om- LDA+DMFT approach): 



PiL) = p{l) 



E 

R 



p\{r)-p'R{r) 



(A6) 



The plane-wave part p(r) can be calculated directly as 
in eq. (|12p from the PAW pseudo wave functions \^ky)- 
For the one-centre terms, the following one-centre density 
matrix may be defined via 



V,. 






(A7) 
Here \pi) are the projector functions from the PAW for- 
malism. With this definition, p\j{r) and p\i{r) can be 
calculated as usual from the partial waves (pi (r) and from 



the pseudo partial waves (f>i{r), i.e. 

pUr.)=J2^^i^*3ir.)Mr). (A9) 

hJ&R 



3. Comparison of the present charge-self-consitent 
scheme to standard LDA+U implementations 

In many implementations of the LDA-I-U scheme, orig- 
inally designed for long-range ordered strongly correlated 
insulators, the local problem is constructed by projecting 
onto a set of angular- momentum channels (i.e., spherical 
or cubic harmonics) within a given range around the cor- 
related site. For instance, the LDA-I-U implementation 
in the MBPP code as well as in the Vienna Ab-initio 
Simulation Package (VASP)^'^ is performed in such a 
way. Here we want to show that the overall interfac- 
ing structure of the present charge-self-consistent scheme 
of extending LDA gives similar results as the traditional 
LDA+U scheme, if a simple purely static mean-field ap- 
proximation to the local interacting problem is used. 

As a simple test, one can compare the variation of 
the total energy resulting from a conventional LDA+U 
calculation with the results obtained from the present 
LDA+DMFT chargc-sclf-consistent calculation now us- 
ing the Hartree-Fock (HF) approximation for the DMFT 
impurity solver. 

We have calculated the equation of state for V2O3 
within the two approaches and the comparison is shown 
inFig.[I3 The LDA+DMFT(HF) (or, simply LDA+HF) 
calculation is performed using the projected local orbitals 
as defined in cq. ([T]) implemented in the VASP code, while 
the results for the LDA+U scheme arc obtained using the 
standard VASP implementation.— The results are rather 
similiar and most of the differences may be due to the al- 
ternative choices for the local projections. Note however 
that especially in the general LDA+DMFT context the 
local projections as defined in cq. ([T]) arc clearly superior 
(e.g. via the resulting wcU-dcfincd local Green's function) 
to the simple angular-momcntum-channcl projections. 
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